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The structural and electronic properties of the wurtzite phase of the InAs and GaAs compounds 
are, for the first time, studied within the framework of Density Functional Theory (DFT). We 
used the full-potential linearized augmented plane wave (LAPW) method and the local density 
approximation (LDA) for exchange and correlation and compared the results to the corresponding 
pseudopotential calculations. 

From the structural optimization of the wurtzite polymorph of InAs we found that the c/a ratio is 
somewhat greater than the ideal one and that the internal parameter u/c has a value slightly smaller 
than the ideal one. In the all-electron approach the wurtzite polymorph has a smaller equilibrium 
volume per InAs pair and a higher binding energy when compared to the zinc-blende phase whereas 
the situation is reversed in the pseudo treatment. The energy differences are, however, smaller than 
the accuracy of standard density-functional codes (~ 30 meV) and a theoretical prediction of the 
relative stability of the two phases cannot be made. 

In order to investigate the possibility of using an LDA calculation as a starting point for many- 
body calculations of excitations properties, we here also present the band-structures of these ma- 
terials. The bands are calculated with and without relativistic effects. In InAs we find that the 
energy gaps of both polymorphs are positive when obtained from a non-relativistic calculation and 
negative otherwise. 

For both semiconductors, we determine the spin-orbit splittings for the zinc-blende and the 
wurtzite phases as well as the crystal-field splittings for the new wurtzite polymorphs. 

PACS numbers: 71.15.Ap, 71.15.Mb, 71.20.-b, 71.55.Eq 



I. INTRODUCTION 

Nanowires (NW:s) are attracting an increasing amount 
of attention in the scientific community because they of- 
fer the possibility of investigating fundamental physical 
properties at the nano scale and because of the numer- 
ous possible applications in electronic and photonic de- 
vices. In particular, InAs-bascd nanowires have already 
been used to improve electronics in ID as, for instance, 
in the case of resonant tunneling diodes [l| and single- 
electron transistors @. These NW:s, typically grown 
on an (lll)B InAs substrate via Chemical Beam Epi- 
taxy (CBE), have a wurtzite crystal structure in con- 
trast to the stable phase at normal pressure and temper- 
ature which is the zinc-blende structure (zb, 3C) with the 
space group F43m (Tj). Despite the successes achieved 
so far with InAs-based NW:s, it has recently been real- 
ized that their wurtzite structure is an important factor 
in the interpretation of experimental data obtained from 
measurements of photo-luminescence (PL) 0, 0, photo- 
current (PC) [f| or other electronic properties 0. 

The situation is similar in the case of GaAs NW:s. 
When the (lll)B GaAs surface is chosen as a substrate, 
the resulting NW:s consist of alternating wurtzite and 
zinc-blend segments Q . Hence, both segments contribute 
to the optical and electronic properties of such a NW. 

Theoretical studies based on total-energy calculations 
[H predict that the wurtzite polymorph (wz, 2H), having 
space group PQ^mc {Cf v ), is a metastable high-pressure 
modification of both the InAs and GaAs compounds. 
Hence, we are faced with new materials, InAs and GaAs 
in the wurtzite phase, for which no previous experimen- 



tal study exists and only little theoretical work has been 
done. 

In this article we investigate the structural and elec- 
tronic properties of these compounds by using all- 
electron calculations based on Density-Functional The- 
ory (DFT) [§]. Indeed, this approach does not require a 
knowledge of material parameters and enables us to com- 
pute numerous properties ranging from the equilibrium 
lattice constant to the cohesive energy of the system. 

DFT in the Local Density Approximation (LDA) 
[ll[ provides an accurate description of the ground 
state electronic structure of solids and such calculations 
are adequate tools for structural studies. The compari- 
son of the calculated band structure to experimental ex- 
citation energies is, strictly speaking, not allowed but, 
nevertheless, provides important information on the elec- 
tronic properties of the material. It is, for instance, well 
known that the experimental band gaps are usually un- 
derestimated by DFT. The excited-state properties of a 
many-electron system, such as the band structure, are 
better described within Many-Body Perturbation The- 
ory (MBPT) as, for instance, by the GW approximation 
[13 , [l3| to the quasiparticle self-energy E. Still, a DFT 
calculation is a necessary step toward the GW calcula- 
tion, since the LDA eigenencrgies and cigenfunctions are 
usually taken to be the zero:th order approximation for 
the perturbative expansion. In particular, the LDA pro- 
vides an accurate Hartree potential which always must be 
treated to infinite order in an extended system. Hence, 
the present band structure calculations allow us to inves- 
tigate how well the LDA can work as a starting point for 
subsequent many-body calculations. 
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Although the LDA cannot describe the excited-state 
properties, it can nevertheless provide a good descrip- 
tion of the properties of the occupied states (3]. This 
means that we can extract useful informations from the 
valence part of the band structure as, for instance, the 
spin-orbit and - for the wurtzite compounds - the crystal 
field splittings. 

Moreover, the all-electron band structures, calculated 
for the first time for the wurtzite phase, offer benchmarks 
for the corresponding pscudopotcntial calculations. 

Related studies reported in the literature are per- 
formed within DFT, use pseudopotentials [H, and 
include the In 4d and the Ga 3d electrons among the 
core states. These approximations lead to an underes- 
timation of the equilibrium lattice parameters and to a 
less accurate description of the band-structure. Here, we 
focus on a study of the ground state properties of InAs 
and GaAs in the wurtzite phase at the DFT/LDA level 
with a full potential description of the atomic species. 
For the calculations of quasiparticlc band structures of 
these compounds we refer the reader to a subsequent pa- 
per \v$. 

The article is organized as follows. After the descrip- 
tion of the theoretical background underlying the calcu- 
lations, we present the results concerning the structure 
optimization of both the zinc-blende and the wurtzite 
phases of InAs. Then, we discuss the band structure 
obtained in the non-relativistic and in the scalar rela- 
tivistic approximations, the latter with and without the 
inclusion of the spin-orbit (S-O) interaction. In the GaAs 
case, instead, we use the experimental lattice constants 
taken from TEM measurements on GaAs NW:s [HI to 
allow for a direct comparison with experiment. Indeed, 
since the deformation potential for the GaAs gap is very 
large [IH, small changes in the lattice constant have a 
strong influence on the band gap. Finally, we present 
the band structures for the GaAs compound calculated 
in the scalar relativistic approximation with and without 
the spin-orbit interaction. 



II. METHODS OF INVESTIGATION 

The electronic structure calculations presented here 
are performed within DFT using the Full-Potential Lin- 
earized Augmented Plane Wave (FP-LAPW) method 1201] 
as implemented in the WIEN2k software package [2lj . 
For the exchange-correlation en erg y in the Local Spin 
Density Approximation (LSDA) El we have cho- 

sen the electron gas data of Ceperley and Alder [22| as 
parametrized by Perdew and Wang [231 ]. 

The APW method [24| is a procedure to solve the 
Kohn-Sham equation based on the division of the unit 
cell into two spatial regions: non-overlapping muffin-tin 
spheres centered at the atomic sites of radius Rmt and 
a remaining interstitial region. The basis set used con- 
sist of atomic partial waves and plane waves in the two 
regions, respectively. Denoting the cut-off wave number 



of the plane waves by Kmax, the number of basis func- 
tions is determined by the product RmtKmax, where 
Rmt is the smallest muffin-tin radius in the unit cell. 
The angular-momentum cut-off of the expansion of the 
wave functions in spherical harmonics within the atomic 
spheres is called Imax- 

Starting from the original APW method, several im- 
proved versions of the method were developed as, for in- 
stance, the Linearized APW method (LAPW) [25|], and, 
with the addition of local orbitals to the basis func- 
tions, the LAPW+LO m, and the APW+lo HI meth- 
ods. Within the WIEN2k code it is possible to use a 
mixed LAPW+LO and APW+lo basis [H| to get the 
advantage of both schemes. For the expansion inside the 
atomic spheres we have chosen the APW+lo basis for 
I = 0,1,2, the LAPW+LO for higher l:s and we have 
used Imax = 10 as the highest / value. 

To ensure that no charge leaks outside the atomic 
spheres we have chosen —9.0 Ry for InAs and —6.0 Ry for 
GaAs as the energy which separates the core and the va- 
lence states. In the InAs case, this choice guarantees that 
the In 4s, As 3d and all higher-energy states are treated 
as band states, i.e., 36 valence electrons per InAs pair. 
Instead, in GaAs, the valence band states start from the 
Ga 3d and the As 3d, leading to 28 valence electrons per 
GaAs pair. 

For both the InAs and the GaAs polymorphs we have 
performed convergence studies to determine the opti- 
mum value of RmtKmax and the number of fc-points 
in the whole Brillouin zone (BZ). We have found that for 
RmtKmax = 11 the total energy is converged to within 
less than 10 meV for both phases. The number of fc- 
points which ensures convergence is 1000 (10 x 10 x 10) 
in the zb and 405 (9 x 9 x 5) in the wz cases, respec- 
tively. The chosen number of fc-points corresponds to 73 
fc-points in the irreducible wedge of the Brillouin zone 
(IBZ) for the zb structure and to 95 fc-points for the wz 
structure. These fc-meshes are used in a modified tetra- 
hedron integration scheme [29[ to determine which states 
are occupied. The Fourier expansion of the electron den- 
sity in the interstitial region {i.e. valence charge density) 
is cut-off at Gmax = 14, leading to 913 and 6524 plane 
waves for the zb and the wz structures, respectively. 



III. STRUCTURE OPTIMIZATION 

The parameters obtained in the convergence study 
were used to carry out the structure optimization of both 
InAs polymorphs in the way described below. 

Due to the cubic symmetry, the equilibrium lattice con- 
stant a z b °f the zinc-blende compound is determined by 
the minimization of the total energy of the system with 
respect to the volume of the unit cell. The energy vs 
volume curve is then fitted to the Murnaghan equation 
of state (30j giving the equilibrium parameters for the 
volume Vq (and hence one lattice constant), the binding 
energy E$ per pair of atoms, the bulk modulus B, and 
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FIG. 1: Deviation from the ideal c/a ratio for InAs in the 
wurtzite phase. 



TABLE I: All-electron structural and energetic parameters 
for the two InAs polymorphs. The equilibrium volume and 
the binding energy refer to one InAs pair. 



wz 
zb 



Voja.u. 3 ) B(GPa) & 



Eo (Ry) 



369.8984 59.6014 4.6170 -16272.35981 
372.2687 61.2977 4.1660 -16272.35959 



its pressure coefficient B' = (dB /dp) p= o . 

On the other hand, the wurtzite phase is character- 
ized by two lattice parameters, namely a wz and c, and 
by the internal parameter u. Consequently, the struc- 
ture optimization of the wurtzite polymorph is carried 
out in four steps ordered according to the importance 
of the parameter under investigation. At first, the equi- 
librium volume is determined by assuming the ideal c/a 
ratio. This gives a first estimate of the equilibrium a wz 
which is then used in the determination of the devia- 
tion of c/a wz from the ideal value. Then, the volume 
optimization is repeated using the new c and a wz finally 
producing the equilibrium values. As the last step, the 
internal parameter u is determined via force minimiza- 
tion. We obtained u/c — 0.3748, which is very close to 
the ideal value (0.375). 

The deviation from the ideal c/a wz ratio is obtained 
by calculating the total energy of the wz polymorph for 
different values of the ratio and then fitting the so ob- 
tained curve to a 4 th order polynomial (Fig. []}. In this 
way, we found c/a = 1.64202, which is slightly larger 
(by 0.5528 %) than the ideal value ^8/3. The resulting 
equilibrium value for c is 7.00118 A. The fact that the 
equilibrium value of c/a is larger than the ideal one sug- 
gests that the wurtzite phase is not the stable poly morp h 
of InAs, in agreement with the rule stated in Rcf. 
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The coefficients of the Murnaghan fit of the energy 
vs volume curves are summarized in Table U for both 
polymorphs. By comparing the results obtained for two 
phases wc find that the equilibrium volume per cation- 
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FIG. 2: Total energy versus volume normalized to one InAs 
pair for the zinc-blende (diamonds) and the wurtzite (cir- 
cles) polymorphs. The continuous lines are the Murnaghan 
fit curves. 



anion pair of the wurtzite structure is smaller than that 
of the zinc-blende by 2.3703 a.u. 3 , corresponding to a 
volume reduction of 0.637 %. 

The result obtained for the binding energy is somewhat 
surprising. Contrary to the observed stability of the zb 
phase, we find that the binding energy of an InAs pair 
is larger in the wz than in the zb case by 3 meV. This is 
clearly seen in the energy vs volume curves, normalized 
to one InAs pair, shown in Figure [2] for both phases. We 
have checked this result by also calculating the total en- 
ergy in the generalized gradient a ppr oximation (GGA) by 
Perdew-Burke-Ernzerhof (PBE) [32| . In the GGA calcu- 
lation we used the same lattice constant and convergence 
parameter as in the case of the LDA. The wz total energy 
was found to be 28 meV lower than that of the zb. We 
note that the difference in the equilibrium total energies 
calculated within the LDA is actually smaller than the 
typical accuracy of calculations carried out within stan- 
dard DFT codes, i.e. ~ 30 meV. We thus conclude that 
the relative stability of the two phases cannot be deter- 
mined within the present accuracy of the calculations. 

As a further check, we have also performed the struc- 
ture optimization by using Ultra Soft (US) pseudopoten- 
tials and plane waves [33|, as implemented in the VASP 
34j. The latter calculations are done by treating 



code 

the In 4rf states as either valence (d-val) or core (d-core) 
states. In the pseudopotential calculations we have used 
a Monkhorst-Pack [35| grid centered on the T point with 
allxllxll and 8x8x8 mesh in reciprocal space for 
the zb and the wz phase, respectively The convergence 
study gives the following kinetic energy cutoffs: 262 eV 
(rf-val, for both polymorphs) and 222 eV and 202 eV (d- 
core) for the zb and wz cases, respectively. The structure 
optimization was performed by minimization of the total 
energy. In the wurtzite case we have assumed the ideal 
values for c/a and u/c. We have thus found that the 
zb phase, having a binding energy lower than the wz by 
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TABLE II: Lattice constants for In As in the zinc-blend and 
wurtzite phases. US denotes Ultra Soft pseudopotentials cal- 
culations with the In 4d electrons included among the valence 
(d- val) or core (d-core) states. The c/a and u/c values were as- 
sumed to be the ideal ones (i.e. 1.633 and 0.375, respectively) 
in pseudopotential calculations. The percentage deviation of 
c/a from the ideal value is 0.5528% in the full potential cal- 
culation. 



search for that value of the internal parameter u for which 
the force on the nuclei vanishes. We searched for this 
point starting from the wz structure with the equilib- 
rium c and a wz lattice constants. The search was carried 
out by means of a minimization routine from the PORT 
library [Z(| choosing 0.5 mRy/Bohr as a tolerance on the 
force. We obtained u/c — 0.37481, i.e. 0.05% smaller 
than the ideal value. 





a zb (A) 


a wz (A) 


C (A) 


c/a 


u/c 


exp. 


6.0542 


4.2839 


6.9954 


1.633 


0.3750 


all-electr. 


6.0428 


4.2638 


7.0018 


1.642 


0.3748 


US d-val 


6.0329 


4.2663 


6.9669 


ideal 


ideal 


US d-core 


5.8023 


4.1060 


6.7051 


ideal 


ideal 



~ 18 meV, is correctly predicted to be the stable phase 
of InAs. We also found that the equilibrium volume per 
InAs pair of the wz phase is smaller than the zb by 0.08%. 

Clearly, the differing results is connected to the differ- 
ent method of calculation, i.e. all-electron versus pseu- 
dopotential calculations. These tests thus provide an 
indication of the accuracy obtainable from calculations 
based on Ultra-Soft paseudopotentials. In our opinion, 
both sets of results must be considered correct in the 
sense that the total energy difference between them is of 
the same magnitude as the overall accuracy in the total 
energies from standard DFT codes. 

The calculated equilibrium lattice constants (in both 
the full and pseudo potential method) are displayed in 
Table [TT| together with the experimental values. The 
experimental values for the wurtzite lattice constants 
are taken from TEM measurements For the zinc- 
blende structure the experimental parameter measured 
at room temperature [361 ] has been extrapolated to zero 
degrees Kelvin. We thus find that the all-electron calcu- 
lations and the pseudopotential calculations with the In 
4d treated as valence electrons give very similar results 
and they slightly underestimate (~ 0.01 — 0.02 A) the ex- 
perimental values for the lattice parameters, as is usually 
the case within the LDA. The results of the pseudopo- 
tential calculation with the In 4c? electrons frozen into 
the core are, instead, rather far from the experimental 
values. This shows the importance of treating the d- 
electrons as valence states in order to obtain a realistic 
description of the InAs compound. In the wz case the re- 
sults reported in the literature [l~5| are calculated using 
HGH pseudopotentials [13] constructed by treating the d 
electrons as core states. The results are: a wz = 4.192 A, 
u = 0.3755, c = 0.6844 A. Also results for the zb case 
are reported in the same article using the same approach, 
leading to a z b = 5.921 A. Other calculations based on the 
all-electron approach ha ve g iven the results a z b = 6.063 A 
[H| and a z b = 6.084 A |39fl in the case of InAs in the zb 
phase. 

The last step in the structural optimization is the 



IV. BAND STRUCTURES 

All the LDA band structures presented in this arti- 
cle are calculated by solving the Kohn-Sham equation 
[h| at the experimental lattice parameters, along lines 
in fc-space which include high symmetry points. We have 
performed such calculations for both the wz and the zb 
polymorphs of InAs and GaAs. 

The motivations for calculating the full-potential LDA 
band structures are as follows. The all-electron bands (i) 
allow us to investigate how well the LDA can function as a 
zero:th order approximation for the perturbation expan- 
sions of MBPT; (ii) allow us to assess the effects of differ- 
ent relativistic approximations; (iii) provide a benchmark 
for pseudopotential calculations based on LDA; (iv) pro- 
vide a good approximation to the electronic structure of 
the occupied states (i.e. valence band states), thus al- 
lowing us to estimate the spin-orbit and the crystal field 
splittings for the new wurtzite structures. 

The WIEN2k code allows for the inclusion of rela- 
tivistic effects at different levels of approximation: non- 
relativistic, scalar relativistic without spin-orbit interac- 
tion and approximate fully-relativistic with the inclusion 
of the spin-orbit interaction. The core states are treated 
fully relativistically 
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unless otherwise specified. In 
the case of the valence states the relativistic effects are in- 
cluded via a scalar relativistic approximation [42| where- 
after the spin-orbit interaction can be added b y us ing the 
scalar relativistic cigenfunctions as a basis In or- 

der to overcome the problem associated with the averag- 
ing procedure inherent in the scalar relativistic approx- 
imation additional local orbitals of definite j character 
have been added to the basis set [45T |. 

When the spin-orbit coupling is taken into account, the 
corresponding Hamiltonian is diagonalized using scalar 
relativistic cigenfunctions in a range of 20 Ry centered 
around the top of the valence band. This is necessary 
because the S-0 splitting in InAs has a strength compa- 
rable to the energy gap and can, thus, not be considered 
as a "small" correction to the scalar relativistic bands. 
In the case of GaAs the S-0 spitting is smaller than the 
energy gap (at low temperature, for the zb phase: 0.33 eV 
and 1.519 cV, respectively [36[), but still the effect is im- 
portant and so we use the same approximation as for the 
InAs. 
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FIG. 3: Full potential band structures of InAs in the wurtzite phase in different relativistic approximations: (a) non-relativistic, 
(b) scalar-relativistic, and (c) scalar-relativistic with the inclusion of spin-orbit coupling. 



A. InAs band structures 

At first, the band structures of InAs in both the zb and 
the wz phases were calculated within the non-relativistic 
approximation. We then found that InAs is correctly 
predicted to be a semiconductor with energy gaps of 
0.395 eV and 0.444 eV in the zb and the wz cases, re- 
spectively. We are grateful to Prof. Almbladh for cor- 
roborating these results using his own LAPW code [46| . 
Hence, within the LDA, the energy gap of the wz poly- 
morph is 49 mcV larger that that of the zb. In the zb case, 
the three p-like states constituting the valence band max- 
imum (VBM) are degenerate at the F point. In the wz 
case two of these states show a crystal field splitting of 
91 meV. The resulting band structure for the wz poly- 
morph is displayed in Fig. [3^. 

The scalar relativistic band structures of the two poly- 
morphs were then calculated without including the spin- 
orbit interaction. Within this approach the InAs energy 
gap is not present anymore. Instead, we found a "nega- 
tive gap" or "wrong band ordering" for both polymorphs. 
Indeed, in the zb case, at the T point, the three p-like 
degenerate states at the top of the valence band have 
a higher energy than the s-like conduction band mini- 
mum (CBM), resulting in an energy gap of T\ c — Ti$ v — 
—0.467 eV. The wz phase exhibit the same problem and 
the calculated energy gap is Ti c — Ti v = —0.359 eV. The 
crystal field splitting can be identified as a splitting of 
84 meV between the the Tq v and T\ v valence band lev- 
els. The scalar relativistic band structure for InAs in the 
wz phase is shown in Fig. [3Jd. 

The main reason for this un-physical result has to be 
ascribed to the notorious limitations of the LDA in re- 
producing band structures. Since the calculated energy 
gaps are typically underestimated, this may result in the 
"negative gap" problem when the semiconductor under 



investigation has a small energy gap. This is indeed the 
case in InAs, for which the experimental energy gap of 
the zb polymorph at zero Kelvin is 0.415 eV [361 ]. 

We should also consider the fact that the binding en- 
ergy of the In Ad states is under-estimated by the LDA 
|47j |. As a consequence the d states are too close to the 
VBM, resulting in a too large repulsion between the p 
and the d states. This has the effect of shifting the VBM 
toward higher energies, hence contributing to the reduc- 
tion of the energy gap. 

Another source of error can be found in the scalar rela- 
tivistic treatment of InAs. Indeed, in this semiconductor, 
the relativistic effects are of the same order of magnitude 
as the energy gap, as can be seen by considering the fact 
that the spin-orbit splitting at zero Kelvin is 0.38 eV [36[. 

We have, therefore, improved our relativistic descrip- 
tion of the system by including the spin-orbit interaction 
as a correction to the scalar relativistic band structure. 
This approach led to spin-orbit splittings in the zb and 
the wz cases of 0.355 eV and 0.336 eV, respectively. In 
addition, we obtained a crystal field splitting of 48 meV 
in the wz polymorph. Unfortunately, the inclusion of 
spin-orbit interaction did not solve the problem with the 
"wrong band ordering" . We obtained the "negative en- 
ergy gaps" —0.453 eV and —0.506 eV for the zb and the 
wz phases, respectively. The band structure with the in- 
clusion of the spin-orbit interaction is shown in Fig. [3ji 
for the wz polymorph. The results for the InAs case are 
collected in Tab. [T[T] 



B. GaAs band structures 

The band structure calculations for GaAs were per- 
formed by usin g th e values of lattice constants mea- 
sured via TEM [18[ on single NW:s having a diameter 
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TABLE III: The experimental and calculated energy band 
gaps, spin-orbit splittings and crystal field splittings (only 
for wz) for the two polymorphs of the InAs. The tabulated 
crystal field splitting is calculated with the inclusion of the S- 
O interaction. All the energies are in eV and the experimental 
values are from Ref. [3t3 |. 





Egap 

non-rel 


Egap 

seal, rel 
w.out S-0 


Egap 

seal rel 
with S-O 


S-O 
splitt. 


Crystal 
field 
splitt. 


zb 
zb (exp) 
wz 


0.395 
0.444 


-0.467 
0.542 
-0.359 


-0.453 
0.415 
-0.506 


0.336 
0.38 
0.355 


none 
none 
48 



of 110 nm. We used the lattice constants a z b = 5.653 A, 
a wz = 3.994 A, c = 6.528 A and u/c = 3/8, i.e. the ideal 
u/c value. As mentioned in the introduction, since the 
deformation potential of GaAs is very strong the calcu- 
lated band gap depends in a critical way on the lattice 
constant. Hence, in order to allow a better comparison 
with experiment, we choose to work with the experimen- 
tal lattice parameters. 

We tested this choice by calculating the relativistic 
band structure (without taking into account the spin- 
orbit interaction) of GaAs in the zb phase with both the 
TEM lattice constant and the equilibrium lattice con- 
stant a z b, eq = 5.6 A taken from Ref. [HI]. The result- 
ing band gaps are 280 meV and 493 meV, respectively. 
Hence, a decrease in the lattice constant of 0.94% results 
in an increase of the band gap of 43.2%. 

The scalar relativistic band structures of both the zb 
and the wz phases of GaAs were calculated using the 
same procedure as in the case of InAs. The main dif- 
ference with respect to InAs is that GaAs is correctly 
predicted to be a semiconductor, even though the LDA 
gap is very small compared to experiment. This behavior 
can be explained with arguments analogous to those used 
in the discussion of the InAs band structures. The band 
structures calculated for the wz polymorph without and 
with the inclusion of the spin-orbit interaction are shown 
in Figs. UandO respectively. 



TABLE IV: The experimental and calculated energy band 
gaps, spin-orbit splittings and crystal field splittings (only 
for wz) for the two polymorphs of the GaAs. The tabulated 
crystal field splitting is calculated with the inclusion of the S- 
O interaction. All the energies are in eV and the experimental 
values are from Ref. [3l| . 
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splitt. 


splitting 


zb 


0.280 


0.167 


0.339 


none 


zb (exp) 


1.629 


1.519 


0.33 


none 


wz 


0.327 


0.211 


0.340 


79 
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FIG. 4: Scalar-relativistic band structure of GaAs in the 
wurtzite phase. 
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FIG. 5: Scalar-relativistic band structure of GaAs in the 
wurtzite phase with the inclusion of spin-orbit coupling. 



Without the S-0 interaction we find that the band 
structure of the wz phase is characterized by a band gap 
of 327 meV and a crystal field splitting of 131 meV. When 
the S-0 interaction is taken into account, the band gap is 
decreased to 167 meV for the zb and to 211 meV for the 
wz phase. The calculated spin-orbit splitting is 339 meV 
and 340 meV in the zb and the wz cases, respectively, 
and the crystal field splitting in the wz band structure is 
79 meV. All these results are collected in Tab. IIVI 



V. SUMMARY AND CONCLUSIONS 

We have reported full-potential all-electron studies of 
InAs and GaAs in both the wurtzite and the zinc-blende 
phases. The structural optimization of both polymorphs 
of InAs is carried out and the result is in good agreement 
with the TEM measured lattice parameters. We have 
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found that, within the present accuracy of standard den- 
sity functional calculations, it is not possible to decide 
which of the phases of InAs is the stable one. Our to- 
tal energy difference between the phases is ~ 3 meV and 
much smaller than the numerical accuracy of the code 
(~ 30 meV). 

In the case of GaAs, the band structures were calcu- 
lated at the experimental lattice constants measured via 
TEM on NW:s having partially wurtzitc structure. 

The LDA band structures have been computed for 
both polymorphs of InAs and GaAs using different rel- 
ativistic approximations. At first, results from the non- 
relativistic approach have been presented for InAs. Then 
we performed calculations using an approximate rela- 
tivistic model in which the core states are treated fully 
relativistically while the valence states are described by a 
scalar rclativistic model. Spin-orbit interaction was then 
added as a correction. When the bands are calculated 
relativistically we find that InAs is wrongly predicted to 
be a zero-gap semiconductor while InAs is a finite gap 
semiconductor within the non-relativistic approximation. 



Instead GaAs is correctly predicted to be a semiconduc- 
tor in the scalar relativistic approximation, even tough 
the calculated band gaps are considerably smaller than 
the experimental one. 

Spin-orbit splittings and crystal field splittings for the 
new polymorph of InAs and GaAs have been calculated 
for the first time. 
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